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SYMBOLS 


isothermal sound speed In collapsing cloud 
visual extinction 
c speed of light 

D distance of observer from accretion disk 

G gravitational constant 

h Planck constant 

blackbody specific intensity 
k Boltzmann's constant 

L shock luminosity 

M mass accretion rate of protostar 

mass of protostellar core 
n hydrogen nucleus number density 

n hydrogen nucleus number density at the shock front 

s 

n^ hydrogen nucleus number density of protostellar wind 

column density of hydrogen nuclei 

Q extinction efficiency of dust 

r distance from protostar 

distance to which grains evaporate 

R distance from protostar to outer radius of accretion disk 
max 

R ^ distance from protostar to surface of cloud 
out 

R distance from protostar to accretion shock 
integrated emission observed from disk 
evaporation temperature of dust 
cloud temperature at collapse 

V speed of accretion shock 

s 

V speed of the wind from a protostellar source 

iii 



A wavelength of radiation 

0 surface mass density of accretion disk 

T optical depth of dust 

Q solid angle of accretion disk as seen by observer 


iv 



SUMMARY 


In this report we estimate the far-infrared and submillimeter continuum and line 
emissions from regions of low-mass star formation and evaluate the capabilities of 
the proposed Large Deployable Reflector (LDR) for studying this problem. The LDR is 
a large space telescope designed for this wavelength range; construction is planned 
for the 1990s. It is concluded that the LDR will be able to probe the temperature, 
density, and chemical structure, and the velocity field in the collapsing envelopes 
of these protostars. In addition, the LDR will be useful in studying the accretion 
shocks on the cores and circumstellar disks of low-mass protostars, and will be able 
to detect shocks driven by protostellar winds. 


1 . INTRODUCTION 


Stairs form through the gravitational collapse of interstellar molecular clouds. 
Although we now possess a wealth of information on the properties of molecular clouds, 
on the circumstances under which stars form, and on the newly formed stars once they 
become visible, little is known about the actual process of star formation. This is 
mainly a result of the lack of an adequate observational tool. As discussed in this 
report, the Large Deployable Reflector (LDR), an orbiting, far-infrared, 20-m-diam 
telescope, will be extremely useful in studies of collapsing molecular clouds. 

In this report we concentrate on application of the LDR to studies of the forma- 
tion of low-mass stars and the possibility of the formation of planets around them. 
First, we review our theoretical knowledge of the star-formation process, emphasizing 
those areas in which the LDR may serve to fill in the large gaps in our understanding 
(sec. 2). Then we describe two types of observational programs possible with the 
LDR: studies of the continuum emission (sec. 3), and studies of the line emission 

(sec. 4) of collapsing molecular clouds. 

The authors wish to acknowledge the support of the Advanced Program Office, NASA 
Headquarters, under RTOP 159-41-01. 


2. THEORY OF STAR FORMATION 


The star-formation process starts with the fragmentation of a parent molecular 
cloud into a number of clumps. The physics of this fragmentation process is poorly 
understood. It is currently believed that turbulence, as well as magnetic fields, 
plays an important role (Myers, 1983; Mestel, 1977). Table 1 lists the average 
properties observed in different types of molecular clouds and clumps in molecular 
clouds. Definite signs of star formation (e.g. , infrared emission from protostars, 

T Tauri stars, maser emission, and compact H II regions) have been detected in all of 
these types of clouds (except for diffuse H I clouds). It should be noted, however, 
that in molecular clouds, the star-formation activity might be confined to the cores. 



TABLE 1.- PROPERTIES OF MOLECULAR CLOUDS 


Type 

L, 

L® 

n, 

cm ^ 

T, 

K 

M, 

M® 

Diffuse H I cloud^ 

5 

20 

80 

2.5x10^ 

Bok globule^ 

.9 

7x10^ 

10 

5.0x10^ 

Small molecular 

5 

8x10^ 

10 

9.8x10** 

cloud‘d 





Dense core‘s 

.1 

>3x10'" 

10 

>3 

Giant molecular 

25 

7x10^ 

10 

1.0x10® 

cloud® 
Warm cores-' 





.3 

>1x10® 

30 

2.6x10^ 


References: “^Spitzer (1978); Martin and 

Barrett (1978); ‘^Snell (1981); ^Myers and 
Benson (1983); ^Solomon (1978); 

/Evans (1978). 


Collapse of a Spherically Symmetric Cloud 

The one-dimensional collapse of a spherically symmetric cloud is relatively 
simple to calculate. As a consequence, much theoretical work has been done on this 
simplified problem, and its characteristics are relatively well understood (see the 
review by Larson, 1978). 

The collapse of a cloud starts isothermically . When the cloud becomes optically 
thick to visible radiation, it heats up. The collapse of a cloud is nonhomo lo gous , 
that is, the center collapses faster than the outside (Larson, 1972). This results 
in a core-envelope structure in which the envelope falls freely in on a core, which 
is already in hydrostatic equilibrium. At the interface, these two zones are sepa- 
rated by an accretion shock. In this shock front, nearly all of the kinetic energy 
of the in-falling material is transformed into internal energy of the gas and is then 
almost completely radiated away. An important difference between the formation of 
low-mass and high-mass stars is that in the formation of high-mass stars nuclear 
reactions commence in the core (e.g., the core reaches the main sequence), while the 
envelope is still accreting. This is not so in the formation of low-mass stars. In 
low-mass stars, most of the generated luminosity is gravitational energy released in 
the accretion shock. It should also be noted that radiation pressure is unimportant 
for the structure of the collapsing envelope around low-mass protostars. 


The velocity of the flow into the accretion shock, Vg , is given by 

1/2 


V 


/2GM /M mil 

= (t7) (^) 


km/sec 


( 1 ) 


The shock luminosity Lg is equal to 


GM M 

L = — ^ 
^S R 


( 2 ) 
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where is the mass in the core, Rg is the radius of the shock, G is the gravi- 

tational constant, and M is the mass accretion rate. For the isothermal collapse of 
a gaseous sphere whose initial state is not far removed from the critical state for 
gravitational stability, the mass accretion rate is approximately given by 

M = J (3) 

where slj is the isothermal sound speed (Shu, 1977; Stabler et al. , 1980). This 
approximation for the mass accretion rate is reasonably accurate until the later 
phases of the collapse when the boundary condition becomes dominant. The shock lumi- 
nosity can now be written as 



(4) 


where Rg = 3x10^^ cm is a typical value for the radius of the accretion shock in the 
collapse of a 1-M® cloud (Stabler et al. , 1980), except for the very first phase 
after the formation of the hydrostatic core. A good approximation of the number den- 
sity n(r) of the freely falling gas is given by 


/R 

n(r) = n^(-f j 


where ng is the density at the shock front. This can also be written as 
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Generally, the collapsing envelope will consist of two zones. In the outer 
zone, the opacity is dominated by dust. In the inner zone, the dust has evaporated. 
The opacity in this zone, caused solely by gas, is four orders of magnitude less than 
in the outer zone. The two zones are separated by a sharp transition region where 
the dust is destroyed. The evaporation temperature of the most-refractory grain 
present in the interstellar medium, graphite, is about 2000 K (Salpeter, 1977). 
Typically this yields for the radius of the transition zone, Rgv> about 10^^ cm. 
Figure 1 shows a sketch of these different zones in a collapsing cloud, and table 2 
lists their approximate sizes. 


The spectral appearance of a collapsing cloud is critically dependent on the 
optical depth, x(A), of the envelope. The principal contribution to the opacity 
comes from the dust: 


T (A) 
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ev 


n^(r) (Tra^)Q(A)dr 
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CLOUD 

BOUNDARY 


Figure 1.- A schematic of the structure of a spherically symmetric collapsing 

cloud. 


TABLE 2.- SIZES OF THE 
DIFFERENT ZONES IN A 
COLLAPSING SPHERICALLY 
SYMMETRIC CLOUD 


Zone 

r, 

cm 

Central core 
Accretion shock 
Dust- free envelope 
Dusty envelope 

3x10^° 

3x10“ 

10^' 

10^^ 


where , ira^ , and Q(X) are the number density, mean surface area, and extinction 
efficiency of the dust, respectively. Using equation (6) and 


we find 


V 1 2, Q(V) 

T n,(r)(^a )- 


N. 


n(r) 


t(X) = 7.6x10' 


/ T R V 
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\10 \l0“ cm/ 
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(9) 


Obviously, a collapsing cloud will be very thick, optically, in the visible 
region of the spectrum during most of the collapse phase. When most of the envelope 
has accreted onto the core, M will decrease, and the protostar will become optically 
visible. (Alternatively, the envelope might also be dispersed by a strong stellar 
wind . ) 
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In contrast to the visual and near- infrared region of the spectrinn, the envelope 
will be optically thin at submillimeter wavelengths during most of the collapse 
phase. Using 


Q(x) ^ if^-3 100 pm 
Q(V) X 


(Whitcomb et al. , 1981) we find 


t(X) 



( 10 ) 


( 11 ) 


It is precisely this point that makes a space-borne, far-infrared and submillimeter 
telescope such an ideal instrument for studying the star-formation process. Figures 2 
and 3 show the results of a numerical calculation of the spherically symmetric col- 
lapse of a l-M© cloud, immediately after the formation of the hydrostatic core, and 
during the last stage of the collapse (Winkler, 1978). From a detailed comparison, 
we see that the approximations discussed above are in good agreement with these 
numerical results. 





Figure 2.- Velocity, temperature, and density distribution of a l-M® protostar 
immediately after the formation of the hydrostatic core (Winkler, 1978). 
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Figure 3.- Velocity, temperature, and density distribution of a I-M 0 protostar at the 
end of the main accretion phase. The hydrostatic core contains 99.3% of the total 
mass (Winkler, 1978). 


Collapse of Rotating Clouds 

Up to now we have only considered the spherically symmetric collapse of a cloud. 
If, however, the parent molecular cloud has angular momentum, then the star-formation 
process will be somewhat different. A collapsing, rotating cloud will form a core 
and a circumstellar disk in the center, surrounded by a freely falling envelope. 
Figure 4 shows a schematic of the different zones in a collapsing rotating cloud; 
table 3 gives their characteristic properties. Note that there will be an accretion 
shock on the disk as well as on the core. As a result, there will be a range of 
shock velocities, depending on the position in the disk. 
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DUSTY 

ENVELOPE 



CORE 


Figure 4.- A schematic of the structure of a collapsing, rotating cloud. 
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TABLE 3.- SIZES OF DIFFER- 
ENT ZONES IN A COLLAPS- 
ING ROTATING CLOUD 


Zone 

R, 


cm 

Central core 

3x10^° 

Dust-free disk 

io“ 

Dusty disk 

io“ 

Dust- free envelope 

10^^ 

Dusty envelope 

10^^ 


The mass in the disk, relative to that in the core, is very sensitive to the 
ratio of the accretion rate to the angular momentum transport rate in the disk 
(Cassen and Summers, 1983). Thus, the thin-disk model of Lin (1981) and the fat-disk 
model of Cameron (1978) for the proto-planetary disk around the early Sun are charac- 
terized by small and large rate ratios, respectively. 

This difference bears directly on the basic problem of planet formation. In a 
thin-disk model, planet formation starts with the condensation of small (submicron- 
sized) dust grains. These dust grains, and possibly interstellar grains which 
escaped destruction, coagulate. The resulting large (centimeter-sized) grains will 
decouple from the gas and settle in the midplane of the disk. The planets are then 
built up from this thin disk of dust (Safronov, 1972). This gradual process of planet 
formation is, however, believed to be slow (~10® yr) , particularly for the outer 
planets. 

In contrast, in a fat-disk model, planet formation is initiated by gravitational 
instabilities in the disk. In this way, dense clumps are formed, which contract upon 
themselves and form the planets (Cameron, 1978). However, this model for planet 
formation also has some drawbacks. In particular, special mechanisms are needed to 
produce nonsolar abundances in the terrestrial planets (Wetherill, 1978). 

Furthermore, it might be difficult to form a rocky core inside the giant planets 
(Cameron, 1984: personal communication). 

Many hydrodynamical studies have been performed on the collapse of rotating 
clouds. Most of these were, however, concerned with studying the development of 
ring, bar, and numerical instabilities. Essentially, no self-consistent numerical 
hydrodynamical calculations have been performed on which the collapse of a rotating 
interstellar cloud is followed all the way to stellar densities. This contrast with 
the one-dimensional case merely reflects the increased complexity of the additional 
dimension. Quite a lot of theoretical effort has also been devoted to the study of 
circumstellar disks around protostars, once the disks have been formed. This work 
is generally based on the viscous accretion disk models of Lynden-Bell and Pringle 
(1974). The surface density and photospheric temperature calculated from such a 
simple analytical model for the solar nebula are shown in figure 5 (Cassen and 
Summers, 1983). For this case, the angular momentum is 5x10^^ g-cm^/sec, a so-called 
"minimum-angular-momentum" nebula (e. g. , the minimum angular momentum found by taking 
the observed angular momentum in the solar system and correcting to solar abundances) . 
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Figure 5.- Surface density and temperature of a "minimum-angular-momentum" solar 

nebula (Cassen and Summers, 1983). The models A, B, and C correspond to times of 
0.5, 1, and 1.6x10® years, respectively. 

Mass-Loss from Low-Mass Protostars 

The presence of strong stellar winds in T Tauri stars was suggested as early as 
1962 on the basis of the P Cygni profile in Ha in some T Tauri stars (Herbig, 1962). 
The best evidence for mass-loss from T Tauri stars comes from molecular observations. 
Carbon-monoxide observations in the vicinity of T Tauri stars and other pre-main- 
sequence, low-mass stars reveal blue-shifted and red-shifted emissions in two opposing 





lobes (Snell and Edwards, 1981). The observed outflow velocities In these low-mass 
stars are of the order of 15 km/sec or less. Emission from shocked, vibratlonally 
excited H 2 has also been detected from the vicinity of low-mass protostars (Simon and 
Joyce, 1983). Evidence for outflow is also seen in the motion of Herbig-Haro (HH) 
objects associated with low-mass protostars. Proper-motion studies and line studies 
of HH objects reveal velocities up to 400 km/sec (Herbig and Jones, 1981; Graham and 
Elias, 1983). The proper-motion vector of the HH objects points back to the exciting 
star. 

The picture that emerges from the wealth of data available, is one of protostars 
losing mass at a high rate. The resultant stellar wind is channeled, quite close to 
the star, into two opposing lobes, possibly by a circumstellar disk. This fast 
(~400-km/sec) wind accelerates the HH objects. The impact on the surrounding molecu- 
lar cloud material causes two shocks to occur. First, the wind will shock against 
the shell of swept-up material (i.e., the wind shock, v ~ 400 km/sec); second, this 
shell of swept-up material will be driven supersonically into the molecular cloud 
(i.e., the molecular shock, v ~ 15 km/sec). 


Large Deployable Reflector 

As discussed before, the main advantage of the LDR in studying star formation 
derives from the fact that the collapsing envelope of a protostar is only optically 
thin for wavelengths longer than about 100 ym (eq. (11)). Currently, the LDR is 
envisioned as a 20-m-diam dish, diffraction-limited to 30 ym. This implies a beam 
size at 100 ym of 1.0 arcsec. At the distance of the nearest site of star formation, 
the Taurus clouds (150 pc), this corresponds to 2.0x10^^ cm (100 AU). Inspection of 
tables 1-3 reveals that the LDR can easily resolve the collapsing envelope. However, 
the core and the disk around it cannot be resolved. Studies of structures in 
proto-planetary disks (planet formation) need a resolution of at least a few astro- 
nomical units. Such a resolution could be achieved with an interferometer with a 
baseline of 1 km. The diameter of each element of such an interferometer could be 
limited to 4 m. A large number of such elements (~25) could then make up for the 
smaller collection area of the individual elements relative to the 20-m design. 
However, linking such an interferometer may present some difficulties. Furthermore, 
such an interferometer will not have the same near-infrared light-bucket possibilities 
as the 20-m design. It should be emphasized that although a 20-m, far-infrared tele- 
scope will not be able to resolve the core and disk, it still will be an extremely 
important Instrument for studies of star formation. In the next two sections, this 
will be illustrated on the basis of possible far-infrared continuum and line emission 
studies with such an instrument. 


3. CONTINUUM EMISSION STUDIES 

Infrared continuum emission is due to warm (T^j ~ 10-10^ K) dust grains. The 
dust in collapsing molecular clouds can be heated to these temperatures through 
either absorption of photons or viscous friction with the gas. The former situation 
arises in the radiative transport in the circumstellar envelope of the protostar, and 
the latter may occur in a viscous accretion disk around the protostellar core. 
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The Envelope 


The spectral appearance of a collapsing nonrotating protostellar cloud with a 
total mass of 3 M 0 is shown in figure 6 . The different models are identified in 
table 4. These spectra have been computed by solving the coupled equations of gas 
and dust hydrodynamics and radiative transfer (Bertout and Yorke 1978). During most 
of the collapse phase of the parental cloud, the spectrum peaks in the far-infrared. 


WAVELENGTH, jum 



WAVELENGTH, pm 



Figure 6 .- The spectral appearance of a collapsing, nonrotating, 3 -M 0 cloud at 
different time steps: models A-G, table 4 (Bertout and Yorke, 1978). 
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Only in the last phase of the collapse (At ^ 10® yr) does the newly formed star become 
optically visible. In effect, t hi s> reflects the high optical depth of the dust in 
the envelope for UV, visible, and near-infrared radiation (sec. 1 and table 4). The 
emission will, therefore, be dominated by the cold dust at the outside, which radiates 
predominantly in the far-infrared. 

The intensity distribution as a function of Impact parameter p at four selected 
wavelengths is given in figure 7 (Yorke, 1980). These are calculated for model F in 
table 4. At a distance of 150 pc (the distance to the nearest region of low-mass-star 
formation in Taurus), the diameter at 30 pm (=2 AU) is about 0.01", which is much less 
than the diffraction limit of the LDR at that wavelength (0.3"). At 125 ym, the 
computed angular size is about 1", comparable to the LDR diffraction limit. 

The calculated total flux at 100 ym is 2><10^^ erg*sec~^ . Assuming an LDR 
detection limit of 2x10“^ Jy at 100 ym at a resolution of v/Av = 10, the LDR will 
make these sources observable out to a distance of about 100 kpc. Besides our own 
galaxy, this Includes the Large Magellanic Clouds (LMC) and Small Magellanic Clouds 
(SMC). At that distance, 1" corresponds to about 0.25 pc, and multiple sources may 
cause some confusion along the line of sight. If it is assumed that the space density 
of T Tauri stars in the LMC and SMC is equal to that observed in the entire Taurus- 
Auriga complex (e.g., about 0.5 pc~^ , Jones and Herbig, 1979) and that the T Tauri 
phase lasts 10 times longer than the collapse phase, a mean separation between col- 
lapsing protostars of about 2.5 pc is obtained. This is a much larger separation 
than the beam size of the LDR and is comparable to the size of dark Lynds clouds 
(~5 pc) in the Taurus-Auriga region. For the denser concentrations of T Tauri stars 
in this complex, the space density is about 8 pc“® (Myers, 1982) leading to a mean 
separation of 1.1 pc. This is still much larger than the diffraction limit of the 
LDR at 100 ym; some confusion along the line of sight may occur in this case however. 


The Disk 

The calculation of the spectral appearance of a collapsing protostar discussed 
above was performed assuming spherical symmetry. In actuality, a circumstellar disk 
will likely evolve (sec. 1). The presence of a circumstellar disk will modify the 


TABLE 4.- CORE AND ENVELOPE PARAMETERS OF A COLLAPSING 
SPHERICALLY SYMMETRIC 3 M PROTOSTELLAR CLOUD AT 
VARIOUS EVOLUTIONARY TIMES AFTER THE FORMATION OF A 
CENTRAL CORE 


Model 

Time after core 
formation, 

10^ yr 

^core ’ 

M® 

L, 

L® 

M, 

10~® Mg/yr 

Ay , 

magnitude 

A 

18 

0.46 

8 

15.6 

4712 

B 

67 

.97 

15 

8.2 

1183 

C 

156 

1.52 

24 

5.0 

580 

D 

332 

2.12 

111 

2.3 

147 

E 

598 

2.50 

90 

.94 

66 

F 

1170 

2.85 

97 

.33 

17 

G 

1345 

2.90 

96 

.20 

.66 
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p, cm 


Figure 7.- Surface brightness as a function of impact parameter for a collapsing, 
nonrotating 3-M0 cloud: model F in table 4 (Yorke, 1980). 


previously discussed results for the radiative transfer, and the disk may be observ- 
able in the far-infrared. We will now calculate the spectrum of a viscous accretion 
disk, neglecting radiative transfer through the circumstellar envelope. This assump- 
tion makes the calculation relatively easy. 


First, it should be noted that the disk is optically thick out to a wavelength 
of about 1 mm. Assuming that the far-infrared properties of the dust in the disk are 
similar to those of dust in molecular clouds (Whitcomb et al. , 1981), we find 



( 12 ) 


for A ^ 30 ym, where o is the surface density of the disk in grams per square 
centimeter. Thus, for a "minimum-angular-momentum" disk (fig. 5) we find an optical 
depth at 1 mm of 4 at a radius in the disk of 10 AU. 


Second, an approximation of the photospheric temperature of the disk is given 
by T = Tg^Rg.^/r, where Tg.^ is the evaporation temperature of the dust .and Rgy is 
the radius in the disk at which the dust will evaporate. The observed Infrared emis- 
sion at a frequency v from such a disk at a distance D is now given by 


S 

V 




I (X)X dX 
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(13) 
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where X is equal to r/Rg^; Rmax outer radius of the disk; and the cosine 

factor takes into account projection effects (Beall et al. , 1984). The specific 
intensity is given by the Planck function at temperature T(X). This yields 


S 

V 





[D(Y ) - D(Y )] 
max ev 


where Y(X) is given by 


and D(Y) is the function 


Y(X) 


hv 

kT(X) 



(14) 


(15) 


(16) 


tabulated in Abramowitz and Stegun (1970). Figure 8 shows the calculated spectrum of 
a proto-planetary disk for Rg^ = 6x10^^ cm (e. g. , the radius of Mercury), 

Umax “ ”^ev " 1^^ ^ (e-g-> the evaporation temperature of silicate material), 

6 = 45° and D = 150 pc (e.g., the distance to Taurus). For comparison, the spec- 
trum of a 4000-K blackbody with a luminosity of 30 L© is also shown. The disk spec- 
trum is much broader than a blackbody owing to the dust temperature gradient. The 
presence of multiple dust components, each with its own temperature gradient, will 
further broaden this spectrum. The total luminosity of the viscous accretion 
disk is about 5 L©. This is, however, highly dependent on the assumed inner radius 
and the condensation temperature. It should be emphasized that the spectrum shown in 
figure 8 will be further modified by radiative transfer in the circumstellar envelope 
As discussed in section 1, the circumstellar envelope will be optically thick for 
wavelengths less than about 100 pm; thus the disk is only observable at far-infrared 
and submillimeter wavelengths. 

During the subsequent evolution of the protostar, the disk has to disappear 
again rapidly (At ~ 10® yr) , because optically visible T Tauri stars do not show the 
excess infrared (10-pm) emiission associated with a viscous disk (Cohen and Witteborn, 
1984: in preparation) . This disappearance can be due to several processes, includ- 

ing angular momentum transport from the protostar to the disk causing the disk to 
expand, accretion of the disk into larger (il cm) objects reducing the surface area 
of solid particles, or disruption of the disk through interaction with a powerful 
stellar wind. Our present knowledge of the star-formation process is insufficient to 
distinguish between these possibilities. It should, perhaps, be noted here that the 
dusty disk around Vega (recently discovered in the far-infrared by IRAS (Aumann 
et al. , 1984) , is probably not a direct remnant of the protostellar accretion disk. 
These particles may be similar to (though possibly larger than) the particles in the 
zodiacal cloud around the Sun. 

Flattened dust distributions have been recently observed in the far-infrared 
looking toward several T Tauri stars, using the Kuiper Airborne Observatory (KAO). 
Figure 9 shows a 100-ym map of the source SSV 13 (Cohen et al. , 1984: in preparation) 
The KAO has resolved this source in one direction. The direction of the outflow of 
the HH objects is perpendicular to the major axis of the flattened dust distribution. 
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Figure 8.- The spectral appearance of a viscous accretion disk (solar nebula). 
For comparison a 4000 K blackbody with a total luminosity of 30 L® is also shown. 


SSV 13 



Figure 9.- A 100-ym map of 
represent 


the source SSV 13, obtained with the KAO. 
Herbig-Haro objects (Cohen et al. , 1984). 


The black dots 
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It is tempting to identify this flattened distribution with the envelope of a col- 
lapsing rotating protostar. With the LDR, it will be possible to resolve such struc- 
tures easily and to study the interaction of a strong stellar wind with the envelope. 


Summary 

Summarizing the discussion on the continuum infrared emission of low-mass proto- 
stars, we conclude that the LDR will be able to detect the IR continuum of a 3-M@ 
collapsing protostar out to a distance of 100 kpc. Although it is projected that the 
Shuttle Infrared Telescope Facility (SIRTF) will be able to detect protostars 
throughout the galaxy, its larger beam size at 100 pm will present some viewing con- 
fusion problems for distances of the order of 10 kpc or larger. Its smaller beam 
size at shorter wavelengths will not be of much help, because the spectrum of a col- 
lapsing protostar peaks so strongly at about 100 pm. For the later phases of the 
collapse (At ^ 10^ yr) , SIRTF will be a good instrinnent for studying the near-infrared 
emission. Although the LDR will not be able to resolve the protostellar core or its 
circumstellar disk, it will be able to resolve the collapsing envelope. 


4. LINE EMISSION STUDIES 


In this section we discuss line emission from the collapsing envelope, the 
accretion shock on the core and the circumstellar disk, and the shocks driven by the 
stellar wind. Each of these is associated with a particular velocity, and high spec- 
tral resolving power (v/Av ~ 10^-10®) is imperative for a meaningful study. 


The Collapsing Envelope 

The gas in the circumstellar envelope of the collapsing protostar is heated by 
gas-grain collisions, and by IR photon excitation followed by collisional deexcita- 
tion. The gas is cooled by emission in molecular rotational and vibrational transi- 
tions. It is important to note that only those regions that are optically thin in 
the dust continuum at a particular wavelength can be probed by emission or absorption 
line studies at that wavelength. 

The far-infrared molecular line emission from molecular clouds during the iso- 
thermal collapse phase has been calculated by several groups (e.g., de Jong et al. , 
1975; Gerola and Glassgold, 1978). Figure 10 shows the calculated infrared spectrum 
of a collapsing cloud with a hydrogen density of 10** cm~^ , a size of 3 pc, and a 
systematic velocity gradient of 1 km* sec~^ • pc"^ for gas temperatures of 20 K 
(fig. 10(a)) and 80 K (fig. 10(b)) (de Jong et al. , 1975; Dalgarno et al. , 1974). 

For comparison, the expected continuum from dust with the same temperature is also 
shown. During this phase of the collapse the cloud is still optically thin in the 
dust continuum, so even the emission lines below 100 pm could be detectable with high 
enough spectral resolution. It is interesting to note that molecules like HCl, HjO, 
and HD can be important coolants of a collapsing molecular cloud despite their low 
abundance (Dalgarno et al. , 1974; Goldsmith and Danger, 1978). 

The intensity of the molecular cooling lines is quite high, which makes them 
detectable even in nearby galaxies. As an example, consider the J = 8 7 transi- 

tion in CO at 325 ym. The computed integrated line intensity in a 3" beam for the 
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Figure 10.- Intensities of molecular rotation lines in a collapsing Interstellar 
cloud at 20 K and 80 K. The line pairs are the different rotational transitions 
of ^'*C0 (left) and ^^CO (right). The other lines are identified in the figure. 
The solid line represents dust emission at the same temperature (de Jong et al. , 
1975). 

80-K case is about 2.5><10“'‘''' W/m^ . That makes this particular line detectable with 
the LDR out to a distance of about 1 Mpc, which includes all galaxies in the local 
group. Here we have assumed a spectral resolving power of 10^, a signal-to-noise 
ratio of 5, and a 25-sec integration time. 

Essentially, no calculations have been performed for the expected cooling line- 
emission of the collapsing protostellar envelope beyond the isothermal phase. Of 
course, the envelope will be optically thick in the visible and near-infrared during 
this phase. Because of the higher gas temperatures in the envelope during this 
phase, an even richer far-infrared line spectrum is expected. The scale size of the 
emitting region is now, however, much smaller cm). This is still larger than 

the diffraction limit of the LDR at 100 ym out to a distance of 1 kpc. It will, 
therefore, be possible to study in detail the density, temperature, and chemical 
structure, and the velocity field in the circumstellar envelope of protostars in the 
Taurus-Auriga clouds during this collapse phase. Such a study will contribute con- 
siderably to our understanding of the physical and chemical processes in collapsing 
molecular clouds. 
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Figure 10.- Concluded. 


The Accretion Shocks 

Besides these thermal emission lines from the circumstellar envelope, emission 
lines from the accretion shock will also be present. In the initial phase, the shock 
velocity is still less than 25 km/sec, and molecules will not be dissociated. Most 
of the kinetic energy of in-fall will be radiated away in the far-infrared rotational 
lines of CO, HD, OH, and HgO. Consider, for example, the collapse of a 1-M@ 
cloud immediately after the formation of the hydrostatic core (fig. 2). The velocity 
of the accretion shock at that time is only 2.5 km/sec. We estimate the total lumi- 
nosity of this shock to be about 2x10“^ L 0 (of eq. (4)). If we assume that this 
energy is radiated by about 100 lines (probably an overestimate), then the flux per 
line from such a cloud in Taurus is about 3x10 W/m^ . With a line width of 
2.5 km/sec, they are much brighter than the continuum. This makes them easily 
observable with the LDR. 

As time increases, the velocity of the accretion shock increases. For shock 
velocities above 50 km/sec, molecules will be dissociated, and the atomic fine- 
structure lines of 01, 0 III, Si II, and Fe II will be the only lines of the shocked 
gas that might be detectable in the far-infrared. However, these lines dominate the 
cooling only below about 5000 K. It is likely that in this phase the temperature of 
the core is higher than this, in which case these lines would be unobservable. Most 
of the energy will be radiated away in optical and UV lines, which are absorbed by 
the circumstellar envelope. 
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Some of the material In the in-falling envelope will rain in on the disk rather 
than on the core. The velocity of the accretion shock on the disk will be much less 
than that on the core. Also, the photospheric temperature of the disk is lower than 
that of the core. Assuming a mass accretion rate on the disk of 10~^ Mg/yr in an 
angular ring with a radius and width of 1 AU, and a shock velocity of 10 km/sec yields 
a shock luminosity integrated over all lines of about 3x10^^ erg/sec. Most of this 
energy is expected to be radiated in the low-lying rotational levels of H 2 O and CO. 
Again, assuming that most of the energy is emitted in 100 lines, that the distance is 
150 pc, and that the line width is 10 km/sec, we find a line intensity of 
10~^ erg* cm“^ *Hz“^ • sec~^ . This is much brighter than the continuum of the disk 
photosphere. Of course, only the far-infrared H 2 O lines will be observable because 
of the high dust opacity in the circumstellar envelope. Again, high resolution is 
needed to separate these lines from the cooling lines of the infalling envelope. 


The Stellar Wind 

Finally, we consider line emission from shocks driven by a stong stellar wind. 

The mass-loss of the protostar will drive two shocks. First, the wind will be shocked 
(v^ ~ 100 km/sec) when it meets the shell of swept-up material; second, because of 
the wind pressure, this shell will expand into the molecular cloud, driving a shock 
of ~10 km/sec. The first shock radiates predominantly in optical and UV lines, 
although some infrared lines will also be present. Table 5 gives the intensity of 
far-infrared lines for a 100-km/sec shock incident on a gas with a density of 
100 cm“^ , and a magnetic field of 1 pG (Shull and McKee, 1979). All of the lines in 
table 4 are detectable with the LDR, as long as the shock fills the beam. Of these 


TABLE 5.- INTENSITIES OF SELECTED 
INFRARED LINES FOR A SHOCK 
VELOCITY OF 100 km/sec, A 
DENSITY OF 100 cm“^ , AND A 
MAGNETIC FIELD STRENGTH OF 1 pG 


Species 

A, pm 

I, 

erg* cm ^ • sec ^ • sr ^ 

C II 

158 

2.1x10“® 

0 I 

63 

1. IxlO"** 

0 I 

145 

1.3x10“® 

0 III 

88 

1.4x10“® 

0 III 

52 

7.4x10"® 

Si II 

35 

4.9x10“® 

S III 

34 

5.5x10“® 


Notes: These intensities are taken 

from Shull and McKee (1979). The 
numerical integrations are termi- 
nated at a temperature of 
500-800 K. Gonsequently , an 
additional 0.06-0.1 eV per H atom 
will be radiated in the fine 
structure lines of Si II, 0 I, and 
C II. 
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lines, the 0 I (63-ym) line is the most interesting. It has been shown that the 
intensity of this line is directly proportional to the mass-loss rate of the 
protostar, 

1(63 ym) = 10"^^ cm~^ sec"^ sr~^ 

where n^ and are the preshock density and velocity of the wind in cgs units 
(Hollenbach, 1984). This makes studies of this line particularly interesting. 

The low velocity of the shell shock ensures that most of the energy is radiated 
in infrared lines. Table 6 gives the intensity of some of these lines for a 10-km/sec 
shock incident on a molecular cloud with a density of 10** cm”^ and a magnetic field 
of 50 yG (Drains et al. , 1983). Again, these lines are bright enough to be observable 
with the LDR, as long as the shock fills the beam. It should be noted that the rela- 
tive intensities of these lines are sensitive to the postshock temperature and den- 
sity and, therefore, to the shock velocity, the preshock density, and the magnetic 
field strength. 

TABLE 6.- INTENSITIES OF SELECTED 

INFRARED LINES FOR A SHOCK VELOCITY 
OF 10 km/sec, A DENSITY OF lO"* cm“^ 

AND A MAGNETIC FIELD STRENGTH OF 
50 yG 


Species 

A , ym 

I, _ 

erg* cm ^ • sec ^ • sr ^ 

01 

63 

2.5x10“'' 

01 

145 

1.2x10“^ 

Cl 

609 

1.7x10"® 

Cl 

370 

1.7x10"® 

H 2 v=0 S(0) 

28 

9.0x10"® 

S(l) 

17 

1 .3x10"" 

S(2) 

12 

3.2x10"® 

S(3) 

9.7 

2.0x10"® 


These intensities are taken from 
Draine et al. (1983). 


One interesting project would be to map HH objects in these lines. The LDR has 
high enough spatial resolution to easily resolve a typical HH object (~10"). High- 
frequency resolution is needed to separate the wind shock from the shell shock, and 
to study the velocity field. In this way, the interaction between the wind and the 
molecular gas can be probed in detail. Among other things, this may help clarify the 
energy and momentum transfer from a stellar wind to a molecular cloud. In particular, 
it may elucidate the nature of HH objects. 


Summary 

We conclude that the LDR will be a powerful tool for studying line emission from 
protostellar objects because of its high angular and frequency resolution. Important 
studies can be made of the temperature, density and chemical structure, as well as 
the velocity field in collapsing envelopes, the accretion shock on their cores and 
circumstellar disks, and of the shocks driven by the stellar wind. 
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